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THEORETICAL INVESTIGATION OF TURBULENT LIQUID-METAL HEAT 
TRANSFER IN CHANNELS WITH HEAT SOURCES IN THE FLUID 

by Robert M. Inman 
Lewis Research Center 

SUMMARY 

A previous analysis to determine the heat-transfer characteristics for turbulent flow 
of a heat- generating liquid metal in a parallel- plate channel (NASA TN D-3473) is ex- 
tended by including the turbulent or eddy conductivity with the molecular conductivity as 
a means of transporting heat. The fully developed velocity distribution is taken as uni- 
form, and an idealized eddy diffusivity function (obtained from Poppendiek) is used. Heat 
is generated uniformly in the flow, and the analysis is carried out for both uniform wall 
temperature and uniform wall heat flux. In solving the energy equation the factor \p, 
which is defined as the average effective value of the ratio of the eddy diffusivity of heat 
transfer to that of momentum transfer, is introduced. The analysis applies in the ther- 
mal entrance region of the channel as well as downstream. 

The effects of the ratios of internal heat generation to wall heat transfer or temper- 
ature driving force on heat-transfer characteristics are investigated. Results for sev- 
eral quantities of engineering interest are presented as functions of the foregoing heat 
generation ratios, a dimensionless axial-distance modulus, and a parameter 

k = fl + 0. 01 xpPrRe®' where Pr is Prandtl number and Re is Reynolds number. 

The solutions for Nusselt numbers are shown to reduce correctly to molecular- 
conduction heat-transfer results as k approaches unity. Results for the fully developed 
Nusselt numbers for liquid metal flow without internal heat sources are compared with 
existing calculations based on the assumption that the eddy diffusivities of heat and mo- 
mentum e jj and are equal and exhibit good agreement (i^ = if/ = 1, where 
* =e H /e M )> The present results apply to systems with Reynolds numbers ranging from 
5000 to approximately 10 6 and Prandtl numbers ranging from 0. 001 to 0. 100. 



INTRODUCTION 


The following analysis is a sequel to an earlier formulation and analytical study of 
the problem of heat transfer to a liquid metal flowing between parallel plates with heat 
sources in the fluid stream (ref. 1). In that study the idealized system which defined 
turbulent heat transfer in a flowing liquid metal was based in part on the following pos- 
tulates: 

(1) The established turbulent velocity profile is represented by a uniform distribu- 
tion (slug flow). 

(2) The thermal eddy diffusivity is small compared to the thermal molecular diffu- 
sivity and is neglected (molecular conduction only). 

The latter postulate implies that the thermal solutions obtained from this model per- 
tain to systems characterized by low Prandtl and Reynolds numbers. These results will 
be referred to in this report as the molecular-conduction heat-transfer solutions. 

The model to be developed in this report improves on the molecular conduction model 
by allowing for turbulent eddying within the liquid metal through the use of an idealized 
eddy diffusivity function (refs. 2 and 3). It is shown in references 2 and 3 that the ideal- 
ized eddy diffusivity function approximates the actual one in the regions nearest the walls 
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over a Reynolds number range of 5x10° to 1x10 . For symmetrical thermal boundary 
conditions, the use of this function should introduce little error, since the heat transfer 
is small in the central regions. For asymmetrical boundary conditions, however, a more 
realistic description may be necessary. 

The assumption is also made in the development of the eddy diffusivity function in 
references 2 and 3 that the eddy diffusivities for momentum and heat are equal. Suffi- 
cient evidence now exists to indicate that, in general, eg and are not equal (refs. 4 
to 10). (Symbols are defined in appendix A.) There is at present, however, no general 
agreement on the form and magnitude of the parameter ip, defined as the ratio 
It is known that the term xp varies across the flow channel as well as with both Prandtl 
and Reynolds numbers. 

To eliminate the radial-dependence drawback, it is convenient (ref. 9) to introduce 
the parameter xp, which is the effective average value of xp across the channel. When 
xp is used, the problem of dealing with local variations in xp is avoided. Reference 9 
considers a method of evaluating the quantity i p for practical use in estimating heat- 
transfer characteristics of turbulent, internal flow of liquid metals. The term \p is 
introduced as a parameter in the present study. No attempt is made, however, to eval- 
uate this quantity beyond assuming xp = 1 when comparing the present results with other 
theoretical predictions wherein the eddy diffusivities for momentum and heat are assumed 
equal. 

The model postulated to approximate the forced convection heat-transfer system 
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under consideration is defined by the following assumptions: 

(1) The established turbulent velocity profile is represented by a uniform distribu- 
tion u(y) = U. 

(2) Longitudinal heat conduction is small compared to convection and transverse heat 
conduction and therefore is neglected. This restriction is shown in reference 11 to in- 
troduce a negligible error for RePr = Pe > 100. 

(3) The eddy diffusivity of momentum distribution varies linearly with distance from 
the channel centerline and as the nine-tenths power of the Reynolds number: 

— = 0. 01 Re 0, 9 ^ 

(4) The internal heat generation is spatially uniform. 

(5) The fluid properties are invariant with temperature. 

(6) A steady state exists. 

The range of validity of this idealized model must be established primarily by com- 
parison with existing numerical and experimental data. The present solutions reduce 
correctly to the molecular- conduction heat- transfer results for Nusselt numbers (ref. 1) 
as the eddy conduction contribution becomes negligible. The fully developed Nusselt num- 
bers in the absence of internal heat generation are calculated on the basis of ip = 1 and 
compared with existing calculations (refs. 12 to 14) and shown to exhibit good agreement. 
Within the knowledge of the author, no heat-transfer measurements for turbulent liquid 
metal flows between parallel plates with heat sources in the fluid stream are available 
for comparison. 

The present results are applicable in the design of liquid- metal- fuel reactors, elec- 
tromagnetic pumps and flowmeters (e. g. , ref. 16), equipment for the electrolytic re- 
moval of gaseous impurities in liquid metals, and liquid metal magnetohydrodynamic 
power generators (e. g. , refs. 17 and 18). 


GENERAL ANALYTICAL CONSIDERATIONS 

A schematic diagram of the system under study is pictured in figure 1, which also 
shows dimensions and coordinates. The direction of fluid flow is from left to right. 

The section of the channel to the right of x = 0 is the region in which the fluid experi- 
ences a uniform internal heat generation Q and to which specific thermal boundary con- 
ditions are applied. Both uniform wall heat flux and uniform wall temperatures are con- 
sidered. The fluid possesses a fully developed velocity profile which is unchanging with 
length. The temperature across the entrance section, x = 0, is taken to be uniform at 
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Figure 1. - Physical model and coordinates. 
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the value t g . For each of the thermal situations, an analysis is carried out to determine 
the heat-transfer characteristics of the fluid at all positions along the duct, that is, in 
the thermal entrance as well as fully developed regions. 


UNIFORM WALL HEAT FLUX 

The starting point of the analysis is the equation expressing conservation of energy 
for turbulent liquid metal flow in a parallel- plate channel (often referred to as a flat duct). 
For fully developed hydrodynamic conditions, this equation may be written as 
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at 

dx 


_d 

9y L v 


(“ +c h) 


at 

3yJ 


_Q_ 

P C r 


( 1 ) 


where a and represent, respectively, the molecular and eddy diffusivities of heat. 
The turbulent velocity distribution is approximated by the uniform distribution 

u(y) s u = constant (2) 


Since equation (1) is linear in the temperature t, it is possible to separate the gen- 
eral problem, which includes both internal heat sources and wall heat transfer, into two 
simpler situations, which are: (1) a problem in which there is a uniform internal heat 
generation Q in the fluid flowing through a duct with insulated walls (q w = 0) (for this 
case, the temperature is designated as tq) and (2) a problem in which there is a uniform 
heat transfer q at the channel walls but no internal heat generation (Q = 0) (in this 
instance the temperature is denoted by t q ). Then, from the linearity of the energy 
equation, the temperature in the combined problem is given by 


t = t Q+ t 


( 3 ) 


The governing equations and boundary conditions for tQ and t may be written as 
follows: 
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and 


u -Q = J- 

9x dy 
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(4a) 
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o 

— — = 0 at y = a (insulated wall) 

sy 


ay 


<-Q > 

— — = 0 at y = 0 (symmetry) 


tQ = 0 at x = 0 (entrance condition) 


(4b) 


at a 

u-2 = J- 

dx dy 


(“- h) 5 ' 


(5a) 


at q 

— 2 = _ at y = a (specified heat flux) 
3y k 


at y- 

— ^ = 0 at y = 0 (symmetry) 
ay 


tq = t g at x = 0 (entrance condition) ^ 


(5b) 


It is convenient to express equations (4) and (5) in terms of dimensionless variables de- 
fined in appendix A; then these equations are given as 


II 

tpl 

a 


9t o| 

ax’ 

3y’ 

\Pr v ) 

ay’_ 


PrK 


(6a) 


at 0 

— — = 0 at y' = 0 and y’ = 1 
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—2 = 0 at y' = 0 
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t = t at x’ = 0 
q e 


(7b) 


It is noted that the dimensionless eddy diffusivity appears together with the reciprocal of 
the Prandtl number in the form of a sum 
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= — + 0. 01 i^Re 0 - 9 (1 - y’) (8) 
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This expression may be inserted into equations (6a) and (7a). The resulting equations 
can be expressed in a simpler form by making the change of variable, 


r = /F 0 + Fj(l - y') (9) 

where Fq = 1/Pr and Fj = 0. 01 xf/ Re 9 ' 9 . The heat-transfer problems to be solved then 
become 


7*Q = l_i , Qa 2 
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— ^ = 0 at r = r n and r n 
9r U i 


tq = 0 at x’ = 0 


(10a) 


(10b) 


6 


J 



(Ha) 



1 1 a 


3x' 4 r 3r 



. at 

1 q 

r 3r 



•N 

at r = r j 
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= 0 at r = r. 


tq = t e at x' = 0 


(lib) 


where the limits are Tq = |/Fq + Fj and r j = |/F^. These equations complete the for- 
mulation of the boundary value problems. The two problems as defined by equations (10) 
and (11) will be analyzed separately. 

The problem of a liquid metal with internal heat sources flowing in an insulated 
channel is considered first. Next, the results are obtained for the problem of uniform 
wall heat transfer with no heat sources. These solutions are then combined linearly to 
yield results for the general case of both internal heat sources and wall heat flux. 


Internal Heat Generation In an Insulated Channel 


A solution for the temperature distribution tq that satisfies, the differential equa- 
tion (10a) may be proposed as follows: 




k Pr 


5s! x+ 


( 12 ) 


where x + = x’/Pr. This form of solution also identically satisfies the entrance and 
boundary conditions given by equation (10b). Equation (12) is a solution of the given 
differential equation and satisfies all the entrance and boundary conditions. It is there- 
fore concluded that the solution for tQ that applies over the entire length of the channel 
is simply 
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( 13 ) 


Qa 2 A 


The local bulk fluid temperature tQ ^(x 4 ) along the channel length, for a uniform heat 
source, is given by 


*Q,b< 0) ] =Qx 


or, since tq b (0) = 0, 


_ Qa- Y + 

b 


(14) 


Uniform Wall Temperature Without Internal Heat Generation 

There is considered next the situation where there is heat transfer q w at the channel 
walls but no internal heat generation. 

For uniform wall heating conditions, it is known that, as the flow proceeds farther 
and farther down the channel, the heat-transfer characteristics asymptotically approach 
a limiting condition. This limit is termed the fully developed heat- transfer situation and 
will be denoted by the subscript d. It is convenient to write the temperature distribution 
t as the sum 




(15) 


Then t* is seen to be an entrance region temperature that is added to t^ ^ to obtain 
temperatures in the region near the entrance of the channel and that approaches zero as 
x becomes large. 


Fully Developed Solution 

Since energy conservation must everywhere be satisfied, it is necessary that t . 

u 

satisfy equation (11a); that is. 
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(16) 
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A characteristic of the fully developed situation for uniform wall heat transfer is that the 
temperature at all points in the cross section rises linearly along the channel. From a 
heat balance on the fluid, it follows that the temperature gradient in the fully developed 
region must be 


9t q, d q w 
9x pUc a 

r 


(17a) 


or alternately 


^ = x+ + G(r) (17b) 

q w a A 

The radial distribution G(r) must, of course, satisfy the conservation- of- energy equa- 
tion (16). When equation (17b) is substituted into equation (16), the governing equation 
for G(r) is found to be 



The boundary conditions on G(r) are found from equations (lib) to be given by 


— = 0 at r = r n 
dr u 


dG 

dr 


> 


— — at r = r 1 

F 1 


(18a) 


(18b) 


J 

The consideration of an overall energy balance on the fluid for the length of channel from 
0 to x. 
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£ [ l q,d< x > - t J dy ' = 


q x 
pUc p a 


leads to an additional condition, upon substitution of equation (17b), on G(r) as 


f * rG(r)dr = 0 
•/r« 


(18c) 


Equation (18a) can be integrated directly. The resulting expression for G(r) is 


C(r)=!!» (lr 2 -rg tor) , 2F 0 (F 0 + F 1 ) 

' 2 0 / „2 


A 


to(Vp 0 + Pi)-l + ^lnh/3 


2F-, 


(19) 


The expression for the fully developed temperature distribution is then given by equa- 
tions (17b) and (19). 


Thermal Entrance Region 

According to equation (15), the temperature distribution in the thermal entrance re- 
gion is found by adding the difference temperature t* to the fully developed temperature 

t j. Since t , has been determined, attention is turned to t*. From the linearity 
q, d q, d q 

of the energy equation, the equation for t* must be 

VI 


5L^1 1 

ax* 



(20a) 


Since the wall heat addition has already been accounted for in the fully developed solution, 

the boundary conditions for t* are given by 

VI 

at* 

— ° = o at r = r 0 and r = r. (20b) 
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and 


‘S-VVd 


or 


> at x 1 = 0 


3 _. 


q w a A 


G(r) 


(20c) 


where equation (17b) has been used in finding the last expression. 
A solution for t* is sought in the form of a product 

4 , 


= 0(r)r(x») 


q w a A 


( 21 ) 


where r and 9 , respectively, depend on x' and r alone. Substituting equation (21) 
into the differential equation (20a) results in two total differential equations. One involves 
t(x’) and the other 0(r). The solution of these equations yields 



q w a A 



+ c 2 Y 0 



( 22 ) 


where /3 is the separation constant arising from the product solution, and Cj and Cg 
are constants in the general solution. 

From the boundary condition 9t*/3r = 0 at r = Tq, the constant Cj is 



Hence the solution for t* can be written as 
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c 2 e 


■/3 2 x* 


J 





( 24 ) 


The constant /3 can be evaluated by substituting equation (24) into the boundary con- 
dition 9t*/3r = 0 at r = r y The solution comprises an eigenvalue problem of the Sturm- 
Liouville type. Solutions fitting this boundary condition can be found only for discrete 
values of j8, that is, /3j, ^ . . . , /3 m , termed eigenvalues. The solutions 
Vg[(2/3j/Fj)r], . . . , V()[(2/ 3 m /F2)r] are the corresponding eigenfunctions. The eigen- 
values are found to be the positive roots of 

J l(w m )Yi(kw m ) - Y 1 (w m )J 1 (kw m ) = 0 (25) 


where w m = (2r 1 /F 1 )0 m and k is a dimensionless diffusivity parameter defined as 


k 


- r 0- J F 0 + F 1 

r l V F 0 



+ 0. 01 \pPrRe 0, 9 


(26) 


The first five values of the roots w m are given in reference 19 for several values of the 

m — 0 9 

parameter k. It is clear that the product i^PrRe * can never be negative, since neg- 
ative values of ^ and/or negative Prandtl and Reynolds numbers are physically impos- 
sible. Therefore, k is equal to or greater than unity for all physically possible values 
of \p, Pr, and Re. It follows from equation (26) that the particular value k = 1 corre- 
sponds to molecular conduction and is interpreted to result for Pr — 0. 

The first five eigenvalues = /^/Iq = w m (k - 1) /4 have been computed from 

the tabulations given in reference 19 for values of k of 1. 25, 1. 67, 2. 50, and 5. 00. An 

2 

analysis of the eigenvalues j3 for k = 1 has been given in reference 1 with the result 
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TABLE I. - LISTING OF EIGENVALUES /3 2 

m 


Diffusivity 


Eigenvalue 


parameter, 

k 

*1 

4 

4 

4 

4 

1. 00 

9.872 

39.48 

88. 83 

158.0 

246. 8 

1.25 

12. 56 

49.99 

112.6 

200.0 

312.4 

1. 67 

17.89 

70. 55 

158.2 

281. 1 

439.4 

2. 50 

32.03 

122.9 

274.2 

485.8 

758.2 

5. 00 

103.3 

373.7 

819.8 

1442 

2242 


= 1) = (mir) 2 ; m = 1, 2, . . . (27) 

The eigenvalues are listed in table I for 
the k values previously mentioned. 

The constants Cg in equation (24) are 
replaced by the coefficients a m corre- 
sponding to the eigenvalues /3 ; hence the 
solution for t* can be written as 


*8 y 

m=l 
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m 
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00 
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(28) 


where 


= 
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J l(kw m ) 


and 




Y l^m) J 0[ w : 


m 



+ 


J lfr"m> Y 0 




(29) 


(30) 


The coefficients A m of equation (28) remain to be determined. When the boundary 
condition at the channel entrance as given by equation (20c) is applied, there is obtained 



(31a) 


and it follows immediately, from the properties of the Sturm- Liouville system for orthog- 
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onal functions, that the coefficients A m are given by 



where G(r) is given by equation (19) and 


V 0 w m 



satisfies the differential equation 


J. _d 
r dr 



(32) 


The analytical evaluation of the coefficients A m has been carried out in detail in appen- 
dix B with the result 


■^m 


4V 0( w m> 

w m^ 2 - [ V 0 (w m) - 


(33) 


It will be useful, in the subsequent analysis, to introduce a new coefficient A m , defined 
as the product A m Vg(w m ): 






w^(k 2 


»> [ y 0< w m> ' 


k^Vj^kw. 


m>] 


(34a) 


TABLE H. - VALUES OF COEFFICIENTS A m 


Diffusivity 


Coefficient 


parameter, 

k 

- A i 

- A 2 

- A 3 

~ A 4 

~ A 5 

1. 00 

0. 20258 

0. 05066 

0. 02251 

0. 01266 

0.00810 

1.25 

. 17870 

. 04494 

. 02000 

.01125 

.00721 

1. 67 

. 14543 

. 03756 

. 01680 

. 00947 

. 00607 

2. 50 

. 09969 

.02771 

.01261 

.00715 

. 00460 

5. 00 

. 04165 

. 01408 

. 00684 

. 00400 

.00261 


The coefficients A m for the particular 
case k = 1 are evaluated in ref- 
erence 1 and given as 

V = k = 1 (Mb) 

(mv) 

Numerical values of A_ are listed 

m 

in table II as a function of the values 
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of k previously mentioned. 


With the / 3 m , V n (w 


m 


, and A m determined, the difference 


is then known from equation (28). It is now possible to proceed with 
tion of the problem. 


temperature, t*, 
the complete solu- 


Complete Solution 


The solution for the temperature t that applies in both the entrance and fully de- 
veloped regions is found by summing the solutions for t^ ^ and t* to yield the result 


= x + + G(r) 

q a//c 
^w 


CO 

I 

m=l 


A m V 0 w m 


-Pm X+ 


(35) 


The local bulk fluid temperature t^ ^ in the presence of uniform wall heat transfer is 
given by 


p uc p^b^-y = v 


or alternately 


q ? - b - - - = x + (36) 

q w a A 


Combined Internal Sources and Wall Heat Transfer 

The solution for the situation where internal heat generation and wall heat transfer 
occur simultaneously can now be written. According to equation (3), the result for this 
combined problem is found by adding the contributions due to each of the separate prob- 
lems to obtain 


t - t e = 


_ ( Qa 2 ^ q w a 


iAV G(r)+ vy 

K K / j 

m=l 


^oKif ) e 


-/3 2 x + 
K m 


(37a) 
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or in dimensionless form, 


t - t 

= (1 + R)x + + G(r) + 

q w a/« 



A m V 0 


m 



where 


R = 


Qa 




(37b) 


(38) 


The heat- flux parameter R is the ratio of internal heat generation to the heat transferred 
at the channel walls and gives a measure of the relative importance, in connection with 
temperature development, of internal heat evolution in the presence of wall heat transfer. 

Of particular practical interest is the wall temperature variation corresponding to 
uniform wall heat transfer and internal heat generation. This quantity can be found from 
equation (37b) by evaluating the equation at r = r^: 




(1 + R)x + + G^j) 



m=l 


(39) 


where G(rj) is given by 



(40a) 


Equation (40a) can be recast to better reveal the sole dependence on the parameter 

k = fl + 0. 01 j^PrRe^' After considerable algebraic manipulation the alternative ex- 
pression for G(r j) is given by 


G(r j) = 


2k 4 In k - k 2 (k 2 - 1) - - (k 2 - lY 
2 

(* 2 - 1) 3 


G(k) 


(40b) 
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A convenient alternate form of equation (39) is obtained by introducing the local bulk 
fluid temperature t^. The bulk temperature for the situation of combined uniform heat 
generation and wall heating is computed from 


or alternately 


t, - t 

^ = (1 + R)x + (41) 

q w a//c 


Then the difference between the wall and bulk temperatures at all stations along the 
channel is given by 



Figure 2. - Wall temperature results for internal heat generation 
with wall heat transfer. 



A noteworthy feature of equation (42) is 
that the variation of the local wall-to- 
bulk-temperature differences is inde- 
pendent of the internal heat generation 
rate. 

The longitudinal variations of the 
dimensionless wall temperatures were 
evaluated from the analytical solution 
(eq. (42)), and the results are plotted 
in figure 2. In interpreting this figure, 
it is important to note that the Reynolds 
and Prandtl numbers appear both in the 
parameter k and in the abscissa. The 
information given in this plot permits 
evaluation of the wall-to-bulk- 
temperature at various stations along 
the channel. 

The fact that the temperature dif- 
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Dimensionless bulk temperature, (t b - t e )/(q w a/K) 


ference t w - t^ (eq. (42)) at any axial location is independent of the heat- flux ratio R 
is believed due primarily to the assumptions of a uniform velocity distribution and con- 
stant physical properties. This is only conjecture, of course, and further investigation 
would be required to support any definitive conclusion. It is worth mentioning, however, 
that the effect of an internal heat source on the fully developed heat-transfer coefficient 
for turbulent flow in a pipe with uniform wall heating has been investigated (ref. 20). 

For flow with a one- seventh power- law velocity distribution and Pr = 0, it is concluded 
in reference 20 that in general the influence of an internal source on the heat transfer 
coefficient (and thus on the temperature difference t w - t^) is slight. It therefore would 
appear that the influence of an internal source on t - would be slight at all axial 
locations, even for Pr > 0, and hence the assumption of slug flow, for example, does 
not introduce a serious error. Again this is only conjecture and should not be considered 
as fact. 

Another quantity which is of practical interest is the fluid bulk temperature variation 
along the length of the channel. The bulk temperature is given relative to the temperature 
of the fluid at the entrance to the channel. The dimensionless bulk temperature given by 
equation (41) has been plotted in figure 3 for parametric values of the heat- flux param- 
eter R. 



Dimensionless longitudinal distance, (x/a)/Pe 


Figure 3. - Bulk temperature results for internal heat generation 
with wall heat transfer. 


Positive and negative values of the 
parameter R are considered in the 
figures. In the present analysis Q is 
taken to be positive (a heat source). A 
positive value of R, therefore, implies 
that q w is positive, or that heat is 
being transferred from the walls to the 
fluid. A negative value of R, on the 
other hand, implies that q w is nega- 
tive, or that heat is being transferred 
from the fluid to the walls. For a posi- 
tive R, therefore, internal heat gener- 
ation and wall heat transfer reinforce 
one another to produce a bulk temper- 
ature larger than that obtained in the 
absence of internal heat generation. 
Conversely, for a negative R, the wall 
heat transfer opposes internal heat gen- 
eration in the bulk temperature develop- 
ment. 

The wall-to-bulk-temperature dif- 
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Figure 4. - Wall temperature ratio in thermal entrance region for 
internal heat generation with wall heat transfer. 


ference for the fully developed situation 
(x — °°) is found from equation (42) to be 

^w “ 

2 = G(k) (43) 

q w a A 

A convenient rephrasing of equation (42) 
may then be carried out by introducing 
the fully developed wall-to-bulk- 
temperature difference, which yields 


V-tb 
(t w - V d 



(44) 


The ratio of temperature differences is plotted in figure 4 as a function of the dimension- 
less axial distance along the channel for parametric values of k by using the numerical 
data listed in tables I and n. Thermal entrance length is commonly defined as the heated 
length required for t w - t^ to approach to within 5 percent of the fully developed value. 

A dashed line is drawn in figure 4 to facilitate finding the entrance length. For example, 
for \}/ = 1, for Re = 56 000, and Pr = 0. 003 (Pe = 168 and k = 1.25), the dimensionless 
entrance length x/a is approximately 8, while for Re = 672 000 and the same Prandtl 
number, Pr = 0.003 (Pe = 2015 and k = 2. 50), the entrance length x/a is approximately 
38. Therefore, the entrance length (for ^ = 1) increases as the Reynolds number in- 
creases. 

It is customary to represent heat-transfer results in terms of a heat-transfer coef- 
ficient h = q w / (t w - t^) and a Nusselt number Nu = hDjj/ k , where Djj is the hydraulic 
diameter (Djj = 4a for the parallel- plate channel). With these definitions, it follows from 
equation (42) that the Nusselt number is given by 


Nu = 


G(k) + 



- -#Lx + 

A e m 
m 


(45) 
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Figure 5. - Nusselt number results for internal heat generation with 
wall heat transfer. 


Numerical values of the Nusselt number 
have been evaluated as a function of the 
dimensionless axial distance x + /4 and 
the diffusivity parameter k and are 
plotted in figure 5. The information 
given in this plot permits rapid evalua- 
tion of the heat-transfer coefficient 
along the length of the duct. 

It is of practical interest to examine 
the Nusselt number in the fully devel- 
oped region Nu- This is obtained 
from equation (45) by considering the 
limit x + — oo : 


Nu d = 


4 

G(k) 


4 (k 2 - if 

2k 4 In k - k 2 (k 2 - 1) - - (k 2 - lV 

2 


(46a) 


Equation (46a) for the limiting case k = 1 assumes the indeterminate form 0/0. Hence, 
it is necessary to apply L' Hospital's rule and to differentiate numerator and denominator 
with respect to k before setting k = 1. Repeated application of this rule yields the final 
result 


Nu d =12 (k = 1) (46b) 

This result agrees with that given in reference 1 for molecular- conduction heat transfer. 
It should be noted that equations (46a) and (46b) are independent of the internal heat gen- 
eration rate. 

The foregoing results are illustrated in figure 6, where the fully developed Nusselt 
number Nu d is plotted as a function of the parameter k. Analytical data from refer- 
ence 13 for fully developed Nusselt numbers employing velocity-profile and eddy dif- 
fusivity of momentum e ^ correlations, and data from reference 14 for fully developed 
Nusselt numbers employing experimentally determined velocity profile and friction factor 
correlations also appear in figure 6. The velocity-profile correlations used in calculat- 
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Figure 6. - Fully developed Nusselt number results for wall heat trans- 
fer and no internal heat sources. Average eddy diffusivity ratio, 1. 


ing eddy diffusivities and Nusselt num- 
bers in reference 13 were based upon 
experimental results of several inves- 
tigators. In addition, the tabulated 
Nusselt numbers (which apply only for 
uniform wall heat flux and no heat 
sources) were calculated in refer- 
ence 13 on the basis of e jj = e ^ and 
were obtained with the aid of a digital 
computer. 

The Reynolds and Prandtl numbers 
to which the Nusselt numbers in refer- 
ence 13 correspond were converted to 
the equivalent values of the parameter 
k (calculated for «// = ip = !)• The re- 
sults are presented for comparison in 
figure 6. The Reynolds numbers 
ranged from approximately 12 500 to 
1. 35x10 with a corresponding range 
in Peclet numbers from approximately 
100 to 10 4 . 


The calculated results presented 
by Martinelli in the discussion section of reference 14 represent corrected values of the 
Nusselt numbers previously presented by Martinelli in reference 12. In these calcula- 
tions, experimentally determined velocity-profile and friction-factor correlations were 
used, as mentioned previously. A value of e jj/e ^ = 1. 00 was also employed in the nu- 
merical computations. The tabulated values listed in reference 14 are for Prandtl num- 
bers ranging from 0 to 1. 0 for Reynolds numbers from 2000 to 10^. However, only values 
in the Prandtl number range 0. 001 to 0. 100 for Reynolds numbers from 2000 to 10 were 
considered here for comparison. 

It can be seen from figure 6 that equation (46a) gives results which are only approx- 
imately 10 percent greater than the fully developed Nusselt values given in reference 13 
over the range of k values shown. The inclusion of a slightly more realistic velocity 
profile in the present study might decrease this difference, but further study is necessary 
before a real conclusion can be drawn. The present analysis for values -of k greater 
than 3. 0 is in good agreement with the numerical data of reference 14. For k less than 
3. 0, the analysis overestimates these numerical data by approximately 10 percent. 

The author knows of neither analytical results including the effect of heat sources 
nor experimental results for heat transfer to liquid metals containing heat sources and 
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flowing turbulently between parallel plates that would allow for quantitative comparison 
with the present model. However, from consideration of the limiting solutions and com- 
parison with existing calculations, it appears that the present relations are reasonably 
correct and can be used to calculate the heat-transfer design characteristics over the en- 
tire length of the passage for the Reynolds number range from 5000 to approximately 10 . 


UNIFORM WALL TEMPERATURE 

Consideration is now given to the problem where a liquid metal at a uniform temper- 
ature t enters a channel whose walls are maintained at a uniform temperature t , dif- 
e w 

ferent from the entrance value (fig. 1, p. 4). In addition, a uniform heat generation 
begins in the fluid at x = 0. The starting point for the uniform wall temperature prob- 
lem is also the energy conservation equation (1), and the velocity profile is again rep- 
resented by equation (2). 

Because of the linearity of the energy equation, it is again possible to write the fluid 
temperature t(x, y) as the sum of two parts, 

t(x, y) = t Q (x, y) + t T (x, y) (47a) 

where t-. corresponds to the problem in which a heat-generating liquid metal at t = 0 
' 0 £ 

enters a channel whose wall temperatures are also t = 0 and t T corresponds to the 
problem in which a nongenerating liquid metal at t 0 enters a channel whose wall tem- 
peratures are t . It will prove convenient to rewrite equation (47a) in the dimensionless 

TV 

form 


t(x, y) - t w _ t Q (x, y) t T (x, y) - t w 


t - t 
e w 


t - t 
e w 


t - t 
e w 


or 


T s Tq + T 


(47b) 


The dimensionless governing equations and boundary conditions for 


Tq and T T are 


”0 F 1 1 8 /_ 8t q\ , Qa 2 
ax’ 4 r 9r y dr J /cPr(t e - t w ) 


(48a) 
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3T Q 

— ^ = 0 at r = r. 
ar 


T q = 0 at r = r j 


Tq = 0 at x' = 0 J 


and 


3Tn 


F 2 
F 1 1 


± (r^\ 


3x' 4 r 3r 

3T m 


0r, 

-\ 


= 0 at r = r r 


3r 


Trp = 0 at r = rj 


> 


T T = 1 at x' = 0 


(48b) 


(49a) 


(49b) 


The solution proceeds in the same manner as for the uniform wall-heat- flux case; 
that is, the two problems are treated separately and the results are then combined to 
yield information for the general situation. For each problem the fully developed situ- 
ation is considered first and then a difference-temperature, applicable in the entrance 
region, is dealt with. These solutions are then superposed. 


Internal Heat Generation With Wall Temperature of Zero 


F ully developed solution . - In the fully developed region the fluid flows isothermally 
(3t Q d /3x = 0) and therefore from equation (48a) tQ ^ satisfies the equation 

OiL (50a 

r 3r \ dr J p 2 «Pr(t e - tj 

with the boundary conditions 
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(50b) 


3T, 




Q,d_ 


3r 


= 0 at r = r. 


Tq, d = 0 at r = r l 


Equation (50a) may be integrated directly. The resulting expression for T Q,d is 

.2 


l Q,d 


Qa 


F 2 K ^e - V 




(51) 


To determine the solution in the thermal entrance region, an entrance temperature 
Tq is introduced such that 


rpsjt HP 

t Q' t Q T Q,d 


(52) 


Since Tq ^ has already been determined attention is now turned to Tq. 

Entrance region solution . - To find the governing equation for Tq, equation (52) is 
introduced into the energy conservation equation (48a). Since T Q,d satisfies equa- 
tion (50a), it follows that Tq must obey 


i 

3x' 4 r 


_ 3 _ 

3r 



(53a) 


The boundary conditions on Tq may be derived from equation (48a), in conjunction with 
equation (52), with the results 


9T n ") 

— — =0 at r = r n 

3r u 

► 

Tq = 0 at r = rj 


(53b) 


and 
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t q = 


-T, 


Q,d 


Qa" 


> at x f = 0 


F 2 KPr ^e " V 


g(r) 


(53c) 


where g(r) = ln(r/r j) - (1/2) (r 2 - r 2 ). 

The solution of equation (53a) that will satisfy equations (53b) and (53c) can be found 
by using a product solution that leads to a separation of variables. This will have the 
form 


rp* _ 2 


Qa' 


2 V ( r\ -*n x ' 


F 2 KPr ^e 
( r \ -2 

|y n — 1 = Vq n and X R are the eigenfunctions and eigenvalues of the new 


(54) 


where V 0 (y„ — 1 = V, 
Sturm- Liouville problem: 


d 2 V ft „ i dV n n 4X 2 

Qi* + ± _°i* + ^V n =0 

2 0,n 


dr 


r dr 


(55a) 


dV, 


V n _ 

= 0 at r = r n ; V n „ = 0 at r = r 1 
dr 0 0’ n 1 


(55b) 


and where y n = (2r ^/Fj)X n . The solution of equation (55a) that satisfies equation (55b) is 


V ° f n *i) " V °’ n = ' Yl(kyn) J ° f n f) + J l (ky n )Y 0 ^n 7^ 


(56) 


where the eigenvalues y n are obtained as the positive roots of 

J l(^n> Y 0<’'n) - = 0 < 57 > 

Values of the roots y n are given in reference 19 for several values of k less than unity. 
Values of k less than unity are physically inadmissible, however, as explained earlier. 
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TABLE HI. - LISTING OF EIGENVALUES 


Diffusivity 

parameter, 

k 

Eigenvalue 

4 

4 

X 2 

X 3 

X 2 

X 4 

X 2 

D 

1. 00 

2.468 

22. 20 

61. 69 

121.0 

200.0 

1.25 

2.851 

27.84 

77. 80 

152.7 

252.7 

1. 67 

3. 555 

38.77 

108.8 

214. 1 

354. 6 

2. 50 

5. 170 

65. 54 

186. 5 

367. 8 

609. 6 

5. 00 

11.48 

186.8 

541.4 

1066 

1784 


It was necessary, therefore, in the present 
investigation to calculate the eigenvalues of 
equation (57) for the values of k cited pre- 
viously. 

2 ”2 

The first five eigenvalues A n = ^/Fq = 

y^(k^ - 1)^/4 are listed in table m for 

values of k of 1. 00, 1.25, 1.67, 2.50, 

2 

and 5. 00. The eigenvalues A n for 
k = 1. 00 are given in reference 1 as 

X^(k= 1) = (2n^l Y. n = 1, 2, . . . (58) 


The next step in the analysis is to determine the coefficients B n of equation (54). 
Applying the boundary condition at the channel entrance as given by equation (53c) results 
in 



(59a) 


According to Sturm- Liouville theory the coefficients B n are then given by 


B n = 



(59b) 


The analytical evaluation of the coefficients B n is carried out in appendix C. The coef- 
ficients B n are obtained from 
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f dV, 


0,n 


4r, 


dr 


r=r. 


B„ = 


n 


/dV, 


0, in 


dr 


r=r-, 


^ky\ 2 

— ) (v 0 n ) 
V r l/ ’ r=r 0 


4r 


I* 

2 n 


n 


(59c) 


It will be useful in the subsequent analysis to define a quantity <p n , given by 




V 1<>V 


_2„2 


- J W*n 


V^(y n ) - k^(ky n ) 1 


(60) 


where 


v o^ y n> = ~ Y i(ky n ) J 0 (ky n ) 


(61a) 


and 


Vi<v„) * - ~~ - Y i<%>-W 

^n ' r / r=r j 


(61b) 


In addition, it is convenient to define a new coefficient B n as 


B n = - 


B n / dV o\ 


dr 


r^(k - 1) \ /r=r j 


4 <P n _ < k ‘ ^n 


(k 2 - l)r 2 


n 


(62) 
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TABLE IV. - VALUES OF COEFFICIENTS B„ 


Numerical values of B are listed in 

n 

table IV as a function of values of the 
parameter k. The values of B n for 
k = 1 are obtained from reference 1 as 

B n = — (k = 1) (63) 

ll O 

*n 


Complete solution . - Now that 
Tq ^ and Tq are known, they can be 
superposed as in equation (52) to obtain the solution that applies over the entire length of 
the channel. This is 


Diffusivity 

Coefficient 

parameter, 






k 

B i 

B 2 

B 3 

1. 00 

0. 81000 

0. 09020 

0. 03242 

1.25 

. 82994 

.08125 

. 02898 

1. 67 

.85528 

.06991 

. 02462 

2. 50 

. 88744 

. 05502 

.01899 

5. 00 

. 90314 

. 03507 

.01150 


^ = 2 r 2 Jjl 

Qa 2 //c(t e - t w ) F x (k 2 - 1) \ r l, 


2 2 
r - r. 




In order to determine, in the later analysis, the local heat-transfer coefficients and 
Nusselt numbers, it is necessary to know the fluid mixed-mean temperature tq which 
may be found as follows: 


A 

T = V _ 

J?CX 


uTqr dr 


ur dr 


r 2 - r 2 Jr 0 
r l r 0 U 


P 


T Q r dr 


since u = U = constant. The temperature distribution from equation (64) is introduced, 
and after integration and rearrangement there results 
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( 66 ) 


T Q, b = 4k 4 In k - (k 4 - 1) - 2(k 2 - l) 2 \ e " X n x+ 

Qa. 2 /K (t e - t w ) 2(k 2 _ l} 3 ~/__j X 2 

n=l 

The situation of a liquid metal containing no heat sources and flowing in a channel with 
/uniform wall temperatures is considered next. This problem is considered in refer- 
ences 2 and 3, but is repeated here for the sake of completeness. 


Uniform Wall Temperature Without Internal Heat Generation 


Fully developed solution . - The temperature Xp(x',r) is the solution to equations (49). 
For large downstream distances, T T ^ approaches zero; that is, t^, ^ — t w . Therefore 
Tp (x' , r) can be obtained from the entrance region solution T^(x', r) defined by 

rp t — m rp 


rp x — rp rp 

a T ~ l T l T, d’ 


Entrance region solution . - The function Trp(x', r) satisfies the equation 


8T T-!i 1± f r <] 


9x» 


4 r Sr 


3r 


(67a) 


with the boundary conditions 


3T* 




= 0 at r = r. 


9r 

T* = 0 at r = r. 


j 


(67b) 


At the channel entrance, the condition is 


T* (0, r) = 1 


(67c) 


The solution to equation (67a) that also satisfies equation (67b) can be written 


OO — Q 

Z /\ -&• 

C n V okfj e 

n=l ' 


( 68 ) 
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where y n = (2r 1 /F 1 )X n and VQ^r/rj) are the same eigenvalues and eigenfunctions ob- 
tained previously from equation (55). The new coefficients C n are found by applying the 
boundary condition at x’ = 0 (eq. (67c)) and using the properties of the Sturm- Liouville 
system with the result 


n 



(69a) 


The integral in the numerator can be evaluated by substituting equation (55) for 
and then integrating. The result of this operation is 



v n r„- 



dr = - 


*j*i 



1 


Hence, when the result from appendix C is used for the integral appearing in the denomi- 
nator of equation (69a), the coefficients C n are given by 


C 


n 



New expansion coefficients C n are introduced for convenience and defined as 


TABLE V. - VALUES OF COEFFICIENTS C fl 


'n 


F 1 C n/ dV 0,n\ 

2 r j \ dr / r=r ^ 


Diffusivity 


Coefficient 


parameter, 

k 

- c i 

- C 2 

i 

Ol 

CO 

- C 4 

“ C 5 

1.00 

2. 0000 

2. 0000 

2. 0000 

2. 0000 

2. 0000 

1. 25 

2. 3661 

2.2617 

2. 2542 

2. 2521 

2.2513 

1. 67 

3.0382 

2. 7021 

2. 6796 

2. 6730 

2. 6708 

2. 50 

4. 5870 

3. 6058 

3. 5388 

3. 5199 

3. 5120 

5.00 

10. 686 

6. 5522 

6. 2278 

6. 1241 

6. 0778 


= (k 2 - 1 )cp n (70) 

For k = 1 the coefficients C n are given 
in reference 1 as 

C n = -2 (k = 1) (71) 

Values of C n are listed in table V for 
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several values of the parameter k. 

The mixed- mean temperature t^ ^ of the fluid is given by its definition 



(72) 


Introducing equation (68) for T T (x', r) into equation (72) and carrying out the integration 
yields 


oo 



n=l 

Combined Internal Heat Sources and Uniform Wall Temperature 


The complete solution for the temperature t(x + , r) that applies for a heat- generating 
liquid metal flowing in a channel with uniform wall temperatures t is found by combin- 
ing the solutions for Tq and T^, in accordance with equation (47b) to obtain the result 


r 


t(x + , r) - t 


t e " *w 


— = T(x + , r) 


Qa^ 


K ^e “ V 






(74) 


When the channel wall temperatures are specified, the longitudinal variation of the 
wall heat flux required to maintain the wall temperatures constant is of practical interest. 
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The heat-transfer rate q at the channel wall according to Fourier's law is given by 


q = 



^1 ^e ~ V /l 3 T n 
2 a \r 3r> 


r=r. 


( 75 ) 


The heat-transfer rate can be found by differentiating the temperature distribution given 
by equation (74) and evaluating the result at r = r j in accordance with equation (75) with 
the result 


qa 


K(t. 


V 



(76) 


where N = QaV#c(t 0 - t ). 

tJ w 

The parameter N indicates the importance of internal heat generation Q relative to 
that of temperature driving force t g - t^. In the present analysis Q is taken to be posi- 
tive, but clearly t - t may be either positive or negative. The parameter N, there- 
fore, may also be either positive or negative. In view of this, illustrative heat-transfer 
results have been evaluated from equation (76) for values of N ranging from -5 to +5 and 
for several parametric values of k. This information is given in figure 7 as a function of 
the longitudinal position coordinate. It is worth additional emphasis that both the position 
coordinate and the parameter k contain the product of Reynolds and Prandtl numbers. It 
can be seen from figure 7 that, for negative values of N, the value of the wall heat flux q 
may change sign at some position along the channel length, which indicates a change in 
direction of the heat transfer at that location along the wall. It is noted, in addition, that, 
for all values of k, the local heat transfer approaches an asymptote that is dependent 
solely on the magnitude of N: 





= N 


This result is obtained from equation (76) by permitting x + to approach infinity therein. 
Thus, for N = 0 (no internal heat sources) the heat flux decreases to zero at large down- 
stream distances. 
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Dimensionless wall heat flux, qa/K(t e - 1 ^) 






.025 


(b) Temperature driving force parameter, -2. 


Dimensionless longitudinal distance, (x/aVPe 


<c) Temperature driving force parameter, 0. 


Figure 8. - Continued. 





(d) Temperature driving force parameter, 2. 


Dimensionless longitudinal distance, (x/a)/Pe 


(e) Temperature driving force parameter, 5. 


Figure 8. - Concluded. 





The fluid mixed-mean temperature for the situation of combined internal heat gener- 
ation and temperature driving force t - t is computed from 

v W 


-r sT b = T Q ,b +T T,b 

w 


where 


oo oo 



H(k) 


4k 4 In k - (k 4 - 1) - 2(k 2 

2 ^ 

2(k Z - 1 ) 


2 

1L 


(78) 


The similarity between equation (77) and the function G(k), given by equation (40b), is 
noteworthy. Then for the fully developed situation (x + — °o) , 




= NH(k) 


(79) 


Local dimensionless mixed-mean temperatures have been evaluated from equa- 
tion (77) and are plotted in figure 8. Also shown are the asymptotic values T^ d as 
given by equation (79). The curves of (t^ - t w )/(t g - t w ) change continuously from an 
initial value of unity to an asymptote of NH(k) as varies from t g to 
t + N(t - t )H(k). In the absence of internal heat generation, the mixed-mean temper- 
ature asymptotically approaches the temperature of the walls. 

A heat-transfer coefficient h and a Nusselt number Nu can now be determined 
from their respective definitions: 


h s 


q 

*b ‘ *w 


hD 


Nu s 


H 


= 4 


qa 


vv 
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^e - V *b " V 



Substitution for qa//c(t e - t^) and (t g - t^)/ (t^ - t ) from equations (76) and (77) leads to 



(80) 


Before equation (80) is considered further, it is of interest to examine the fully developed 
heat-transfer conditions. For nonzero values of N, the fully developed condition is ob- 
tained from equation (80) when x + — °° and the Nusselt number Nu^ is given by 


N u d 


_4_ 

H(k) 


(81) 


It is noteworthy that this result is independent of the parameter N. For molecular- 
conduction heat transfer, k = 1, H(l) = 1/3, and the fully developed Nusselt number is 

Nu d =12 (k = 1) (82) 


which agrees with the result given in reference 1. 

The Nusselt number at any axial location for N = 0 is obtained from equation (80) as 
the simpler expression 



(83) 


The fully developed heat-transfer condition in the present situation occurs for x + suf- 
ficiently large- so that only the first term of each of the series need be considered. The 
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60i 


Diffusivity 

parameter, 



10 1 I I I l i 1 i i 

0 . 025 . 050 . 075 .100 

Dimensionless longitudinal distance, (x/a)/Pe 

(e) Temperature driving force parameter, 5. 

Figure 9. - Concluded. 

fully developed Nusselt number then follows as the very simple expression 

Nu d = (84) 

For k = 1 (molecular- conduction heat transfer), 

Nu d (k = X) = 4(|^ 2 

= 9.87 (85) 

in agreement with the result given in reference 1. 

Figure 9 shows the variation of the local Nusselt number as a function of the longi- 
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tudinal position coordinate for values 
of N ranging from -5 to +5, with k 
appearing as the family parameter. 

An overall inspection of the figure re- 
veals that, for positive values of N 
(including the case where N = 0), the 
local Nusselt numbers are large near 
the channel entrance and decrease con- 
tinuously with increasing downstream 
distance. The Nusselt numbers plotted 
in figure 9 for negative values of N, 
however, indicate that the definition of 
the heat- transfer coefficient can become 
meaningless in these instances. This is 
due to the fact that the mixed-mean fluid 
temperature can become zero, or can 
have a sign different from that of the 
wall heat flux. Thus, zeroes and infin- 
ities can exist in the Nusselt numbers. 

Numerical values of the fully devel- 
oped Nusselt number have been evalu- 
ated from equations (81) and (84) as a 
function of the parameter k, and are plotted in figure 10. Inspection of the figure re- 
veals that the effect of internal heat generation is to raise the Nusselt number above the 
value attained by a nongenerating liquid metal. It is interesting to compare figure 10 
with its counterpart for uniform heat flux, figure 6 (p. 21). These figures portray sim- 
ilar trends and, in addition, the numerical values are not drastically different. 

The fully established Nusselt number results, in the absence of internal heat sources 

(eq. (84)), have been compared to the values calculated by Martinelli for uniform wall 

heat flux and presented in the discussion section of reference 14. The tabulated 

Nusselt number values in the Prandtl number range of 0. 001 to 0. 100, for Reynolds num- 
4 6 

bers from 10 to 10 , were considered here for comparison. The comparison (fig. 11) 
illustrates that the two sets of calculations are in fairly good agreement, with the 
uniform-wall-temperature Nusselt numbers falling about 10 percent below the uniform- 
wall-heat- flux Nusselt numbers. 

Within the author's knowledge there are no uniform, symmetrical wall-temperature 
analyses for liquid metals flowing turbulently between parallel plates for comparison. 
Reference 15 considers heat transfer to a turbulently flowing fluid between parallel walls 
with asymmetric, uniform wall temperatures. It is concluded in reference 15 that, for 



Figure 10. - Fully developed Nusselt number results for uniform wall 
temperature without internal heat sources and for uniform wall 
temperature with internal heat generation. 
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Figure 11. - Comparison of analysis for fully developed Nusseit num- 
ber for uniform wall temperature without internal heat sources to 
numerical data obtained from reference 14. 


the case of equal wall temperatures, the 
solution given therein can be shown to 
coincide with the solutions given in ref- 
erence 12 for the case of equal wall tem- 
peratures; Reference 12, however, 
treats the situation of heat transfer to 
fluids flowing between parallel plates 
with uniform wall heat fluxes and the 
conclusion drawn is therefore not clear. 

Presently available experimental 
data are not sufficiently extensive to 
establish the accuracy of the Nusseit num- 
bers given by equations (84) and (85). 

Within the author's knowledge there 
are also no heat-transfer analyses for 
heat- generating liquid metals flowing 
between parallel plates that apply over 
the entire length of the duct and that con- 
sider simultaneous transfer of heat by 
molecular conduction and eddy diffusion. 
There are, in addition, presently no ex- 
perimental results for liquid metal flow 


between parallel plates with internal heat generation with which to compare the theoret- 


ical results. 


CONCLUDING REMARKS 

Solutions have been obtained for heat transfer to liquid metals flowing turbulently 
between parallel plates with heat sources in the fluid stream. In order to provide math- 
ematical simplicity, the slug- flow model for the turbulent velocity distribution was uti- 
lized for the analyses. The effect of transverse eddy diffusivity variations, of importance 
for high Reynolds number systems, has been included through the use of an idealized eddy 
diffusivity function. This function approximates the actual one in the heat-transfer re- 
gions of importance. 

The present analyses allow for consideration of Reynolds numbers ranging from 
5xl0 3 to approximately 10®, in virtue of consideration of transverse heat flow by turbulent 
eddying as well as by molecular conduction. Various quantities of engineering interest for 
both uniform wall heat flux and uniform wall temperature cases are determined in both 
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the thermal entrance and fully developed regions as functions of ratios of heat generation 
rate to wall heat transfer or temperature driving force, a dimensionless axial distance 
modulus, and a diffusivity parameter k. It has been shown that the effects of Reynolds 
and Prandtl numbers can be expressed concisely by the diffusivity parameter. 

Evidence presented indicates that the fully developed Nusselt number solutions re- 
duce correctly to known molecular- conduction heat-transfer results as k approaches 
unity. The present solutions, moreover, are shown by comparison to give satisfactory 
accuracy in the limiting case of no internal heat generation in the liquid metal. A truly 
quantitative comparison would require, on the theoretical side, inclusion of slightly more 
accurate velocity and eddy diffusivity descriptions in the boundary value problems and, 
on the experimental side, experimental data for the case in which there is internal heat 
generation in the flowing liquid metal. 

The results indicate the effects of turbulent or eddy diffusivity and internal heat 
sources on the complicated problem of heat transfer to a flowing liquid metal, and it is 
believed that the simple engineering approach to the problem demonstrated here should 
be adequate for preliminary design purposes. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 20, 1966, 

129-01-09-07-22. 
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APPENDIX A 


SYMBOLS 


V 

coefficients in series expansion of 
temperature for case of wall 

F 1 


heat transfer and no internal 
a 

heat sources, — — - 

J l< kw m> 

G(r) 

V 

coefficients defined by product 



G(k) 

a 

half- height of channel 

g(r) 

a m 

coefficients in series expansion of 

H(k) 

temperature for case of wall 
heat transfer and no internal 

h 


heat generation 

I< 

B n 

coefficients in series expansion of 

1 


temperature for case of internal 

h 


heat sources and zero wall tem- 

h 


perature 

h 

*n 

coefficient defined by eq. (62) 

J o 

C n 

coefficient in series expansion of 


temperature for case of uniform 
wall temperature and no internal 
heat generation 

J i 

- 

k 

C n 

coefficient defined by eq. (70) 


C 1 

constant of integration 


c 2 

constant of integration 

N 

C P 

specific heat of fluid at constant 
pressure 

Nu 

d h 

hydraulic diameter of duct, 4a 

F 0 

reciprocal of Prandtl modulus, 

Pe 


1/Pr 

Pr 
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function of average value of x}/ and 

— 0 9 

of Reynolds modulus, 0. 01 i/'Re ’ 

transverse temperature distribution 
in fully developed region for case 
of wall heat transfer and no inter- 
nal heat generation 

function of k given by eq. (40b) 

function of r defined by eq. (53c) 

function of k defined by eq. (78) 

local heat-transfer coefficient, 

- V or V'fV - *b> 

integral defined by eq. (B7) 
integral defined by eq. (B12) 
integral defined by eq. (Cl) 
integral defined by eq. (C6) 

zero order Bessel functions of 
first kind 

first order Bessel functions of 
first kind 

diffusivity parameter, 

r Q /r 1 = Vl + 0. 01 ^PrRe 0, 9 

ratio of internal heat generation 
rate to temperature driving 
force, Qa 2 //c(t e - t w ) 

local Nusselt modulus, hD|j//c 

Peclet number, RePr 
Prandtl number, v/a 


mi ■ ■■■■■■■I ■ i i n mi n ■■■iiiii ■ him 1 1 ■ i 



Q 

heat generation rate/volume 

q 

local heat-transfer rate per unit 
area at channel wall 


prescribed wall heat flux 

R 

ratio of heat generation rate to 
wall heat flux, Qa/q w 

Re 

Reynolds number, 4aU/^ 

r 

transformed coordinate, 
+ Fj(l - y’) 

r 0 

value of r at y* = 0, ^Fq + Fj 

r l 

value of r at y* = 1, ^Fq 

T 

dimensionless temperature, 
(t - t w )/(t e - t w ) 

t 

fluid temperature 

U 

uniform velocity in x-direction 

u(y) 

fluid velocity in x-direction 

v o 

linear combination of zero order 
Bessel functions, eqs. (30) and 
(56) 

V 0, n 

linear combination of zero order 
Bessel functions, eq. (56) 

V 1 

linear combination of first order 
Bessel functions, eq. (61b) 

w m 

modified eigenvalue, ( 2r i/Fj)0 m 

X 

longitudinal coordinate measured 
from channel entrance 

x f 

dimensionless longitudinal coordi- 
nate, 4x/aRe 

x + 

dimensionless longitudinal coordi- 
nate, x’/Pr = 4x/aPe 

Y o 

zero order Bessel functions of 
second kind 


Yj first order Bessel functions of 

second kind 

y transverse coordinate measured 

from channel centerline 

y' dimensionless transverse coordi- 

nate, y/a 

a molecular diffusivity for heat, 

K /P° p 

0 separation constant 

0 eigenvalues defined by 

0 m = (F 1 /2r 1 )w m 

^m eigenvalues defined by 0^ = /3 ^/Fq 

y separation constant 

y n eigenvalues of eq. (57), (2rj/F])X n 

e eddy diffusivity 

e jj eddy diffusivity of heat 

e jyj eddy diffusivity of momentum 

0 function of r, eq. (21) 

k thermal conductivity 

X eigenvalue 

X n eigenvalues defined by 

\=(F 1 /2r 1 )y n 

rt 9 

X^ eigenvalues defined by X o /Fq 

v kinematic viscosity 

p fluid density 

t function of x, eq. (21) 

cp n function defined by eq. (60) 

ip ratio of eddy diffusivity for heat 
transfer to that for momentum 
transfer, 

ip average value of xp 
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j// n function defined by eq. (59c) 
Subscripts: 

b bulk or mixed-mean condition 

d developed 

e entrance 

Q internal heat generation, insulated 

wall or zero wall temperature 


q wall heat flux, no heat generation 

T prescribed wall temperature, no 

internal heat generation 

w wall 

Superscript: 

* entrance region 
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APPENDIX B 


EVALUATION OF SERIES COEFFICIENTS A 


m 


The coefficients A m of the series expansion (eq. (28)) are determined as the quo- 


tient of two integrals (eq. (31b)): 


X lrG<r)v °( wm ^) 

f lrvi ^y 


(31b) 


In order to be able to solve for the coefficients A m as the quotient of the two known in- 
tegrals, it is necessary that the sequence of functions Vq ( w m — j form a system orthog- 

V r l / 

onal with respect to the weight function r over the interval (rg, r j). This means that 
the functions have the property 


/' 


rV 0( w m — ) V o( w n— W = 0 m * n 


(Bl) 


The orthogonality of the set of characteristic functions Vq (w m — j will now be demon- 

V r V 


strated 


Let V 0 |w m — J s V 0 m and V 0 |w n — j s V Q n be the solutions associated with 


two distinct values of / 3 , namely, and / 3 ^. This means that 

d2 


dV r 


4# 


m 


1A rAAll + 

rdr \ dr / p2 


V 0,m = ° 


(B2) 


with dV Q m /dr = 0 at r = r Q and r = r^, and 
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1 J_ 
r dr 



(B3) 


with dV Q n /dr = 0 at r = r Q and r = rj. Equation (B2) is multiplied by V Q ^ n and 
equation (B3) by V 0 m , and then equation (B3) is subtracted from (B2). The result, 
after transposing, is 




rV n V n 
0, m U, 





dr 


- V, 


0, m 


dV 0.n\ 

* / 


If this equation is integrated between Tq and r^, there is obtained the result 

r l 

(B4) 

r 0 

The right side of equation (B4) is zero in virtue of the boundary conditions cited previ- 
ously, and therefore 




rVn ™ v n n dr = 0 
0, m 0 ? n 


(B5) 


— — — — o “2 

Since j8 and 0 n are two distinct values of /3, the difference /3 m - /3 r cannot vanish. 

Hence 



0 n * m 


(B6) 


and orthogonality of the functions is demonstrated. 

Consider again equation (31b). Let Ij equal the integral appearing in the denom- 
inator: 
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(B7) 



dr 


The functions 



satisfy the differential Bessel equation 


_d 

dr 



(B8) 


If this equation is multiplied by 2r(dV Q /dr), and the resulting equation is integrated by 
parts over the range r q to r^, there results 



Into this expression there are substituted the boundary conditions 


which give 


dV Q 

— — = 0 at r = r n 
dr u 


and r j 



(BIO) 


and thus 
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"1=7 [ V 0 (Wm> ' k2v O (kW m>l 


(Bll) 


Consider next the integral appearing in the numerator of equation (31b): 




rG(r)V 0 (w m — Jdr 


(B12) 


The integral I 2 can be evaluated by substituting G(r) from equation (19), integrating by 
parts, and utilizing equations (B8) and (BIO). This process gives finally 


J 2 = 


2/3 


2 L 
n 


(' 2 - 4 ) 


Vr 


2r l V 0 ( w m> 

/ 2 2\ 2 
( r l ' r 0/ m 


(B13) 


Therefore, when equations (31b), (Bll), and (B13) are used, 


A„ = 

m 


4r l V 0( w m) 


( 2 2 \ 2 

( r l - r 0/ w m 


v §< w m> - k2v ?0™m) 


Substituting k = r^r ^ yields 


Am = 


w 2 (k 2 - 1) 
m 7 


4V o( w m> 

Vjj<w m > ’ k2v 0< kw m> 


(B14) 
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APPENDIX C 


EVALUATION OF SERIES COEFFICIENTS B 


n 


The coefficients B n of the series expansion (eq. (54)) are determined as the quo- 
tient of two integrals (eq. (59b)) : 


B n = - 


/ 


rg(r)V Q (y n — ]dr 


£-^4 


(59b) 


Idr 


As was the case for the sequence of functions VQ^w m — j, the present functions 

V 0 (y — ] form a system that is orthogonal with respect to the weight function r over 
\ r l / 

the interval (rg, r j). The proof follows the procedure given in appendix B and need not 
be repeated here. Hence, it is possible to solve for the B R in the manner prescribed 
by equation (59b). 

The general solution to the integral 


■£•44 


dr 


(Cl) 


may be obtained by the technique used in solving the integral Ij in appendix B. The re- 
sult is 




, 4 v 2 

* 2 1 


frVo ) 2 


r l 


r 0 


8F 0 X 2 


*1 


(C2) 


53 


Into this expression there are substituted the boundary conditions 


dV c 

dr 


= 0 at r = r Q ; Vq = 0 at r = rj 


This substitution yields 



4F Q X 2 r2 



(C 3) 


(C4) 


and therefore 



where g(r) = Tq lnfr/r^ - (1/2) (r 2 - 
g(r) from this expression and Vq fy n 



The integral can be evaluated by substituting 
from equation (55a), namely, 




4 Vn 


r dr 



integrating by parts, and utilizing equation (55b). There is obtained the result 
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I 

f 

f 


*4“ 


4 Vn 




fdV 


0, n 


dr 


r=r . 


-■o< v o,„> r 

? r=r f 


F i r i 

2 vS 


/dV, 


0,n] 


dr 


+ r 0< V 0,n> _ 


r=r. 


r=r. 


(C7) 


Since g(r^) = 0, however, equation (C7) may be rewritten in the simpler form 


>4 = 


F t r l / dV 0,n\ 


8F ^n 


dr 


r=r. 


k iV, 


0,n| 


8 ViV dr 


(C8) 


r=r. 


Therefore, using equations (59b), (C5), and (C8) yields 



(C9) 
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